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Abstract 

In the present letter we employ the method of the dynamical renormalisation group 
to compute the average mass distribution of aggregating point particles in 2 dimensions 
in the regime when the effects of local mass distribution fluctuations are essential. 
PACS numbers: 05.10.Cc, 82.20.Mj. 
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The statistical properties of systems of aggregating diffusive particles have been studied 
using methods of kinetic theory since the pioneering work of Smoluchowski, [Q. The kinetic 
equation derived in this work almost a hundred years ago has proved extremely useful in 
describing properties of various aggregation systems, see Q for a review. However, as it 
became clear in the 80's, the Smoluchowski equation provides a sort of mean field theory 
description of the process of aggregation and is therefore valid if only the fluctuations of 
particle density around mean values are small. 

Both theoretical analysis within the mean field theory Q and numerical experiments 

discovered the existence of the critical dimension dc of the ambient space which marks 
the transition from non-fluctuating large time statistics of a given model of aggregation to a 
fluctuation-dominated one. This means in particular, that the scaling properties of the large 
time asymptotics of average mass distribution in d > dc are well described by an appropriate 
Smoluchowski equation, whereas in d < dc the mean field theory answer for the large time 
asymptotics of the average mass distribution becomes incorrect. 

By analogy with equilibrium critical phenomena one should expect that mean field theory 
answers receive logarithmic corrections in d = dc- These have indeed been observed in 
numerical experiments Q and explained heuristically Q . 

However, the theory of stochastic aggregation which would enable one to account for 
fluctuation effects systematically, in particular to derive the above mentioned logarithmic 
corrections in d = dc, is still absent. 

The aim of the present letter is to build such a theory for the simplest model of stochastic 
aggregation of diffusing point particles characterized by constant diffusion and reaction rates. 
This model was originally introduced in Q . Its critical dimension was found to be equal to 
2. 

We start the discussion by defining the model in terms of the master equation. We then 
use Doi's formalism to derive the functional integral representation for the average mass 
distribution. We then establish the connection between the resulting effective field theory 
and the well understood effective field theory of A -I- A ^ ^4 model pi . The connection is 
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used to compute the average density in 2 — e dimensions. This allows an easy derivation of 
logarithmic corrections to the mean field theory answers in d = 2, the determination of the 
exact scaling properties of the answer in d = 1 and a discussion of the mean field theory 
applicability in d > 2. 

Consider a set of point particles on a d-dimensional lattice Z'^ each performing a random 
walk characterized by a rate D. 

Suppose that each particle carries an integer charge q. Any two particles with charges qi 
and q2 positioned at a given lattice site can coalesce with probability A per second to create 
a particle with a charge qi + q2 in the same position. 

It is convenient to think of our particle system as living on a (l + d)-dimensional lattice T 
with d space coordinates and one charge coordinate. We assume that initially the occupation 
numbers of particles of a given charge at different lattice sites are independent identically 
distributed Poisson random variables. We also assume that at t = there was no particles 
with negative charges. Then there will be no such particles at any later time and our problem 
becomes identical to the problem of mass aggregation. 

We wish to study the time evolution of the average number of particles at a site of F 
which corresponds, after taking the continuum limit, to the mass distribution of particles. 

A microstate of our particle system is defined by specifying the number of particles (the 
occupation number) at each site of F. The evolution of the probability measure P on the 
space of microstates is governed by the so called Master Equation. Our further considera- 
tions are based on the following two facts about the Master Equation: 

(i) The Master Equation is hnear; 

(ii) The Master Equation is the first order differential equation w. r. t. time. 

Thus following the already standard construction of |^ one can construct a path integral 
representation of the solution to the Master Equation and, consequently, of any correlation 
function. Leaving details of the derivation to j7| we present here the final path integral 
expression of the average mass distribution in the continuum limit and study it. 

Let Ax be the lattice spacing in the coordinate space, Aq - in the mass space. The 
continuum limit is taken according to the following rules: 

^{Axf ^ Xo, DiAxf ^ D. (1) 

Then the average density is 

n{q, t)^ j D^{-K, q, t)D(^* (x, 9, t),/)(x, q, t)e-^l'^''^'^+I 1 d<ino{<iW (x,,) ^ 

where 

S = dtl^J d'^xdq(f>*{x,q,t)dtcb{^,q,t) + H[(j,,<l>*]j (3) 
is an effective action, 

H = D J d'^x J dqV(t)*{x, q) ■ V(/)(x, q) - Xq J d'^x J J dqdq' 

(^0*(x, q + q')- 20* (x, q) - 0*(x, g)0*(x, g')) 0(x, 9)0(x, q') (4) 
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is an effective Hamiltonian; no(q) is the initial average mass distribution, 



^^(9,0 = ^96 « (5) 

if g > and otherwise. Here uq is the total initial particle density, M- the total mass 
density. 

The non-hermiticity of the effective Hamiltonian reflects the irreversibility of the process 
of aggregation. As in case of the equilibrium statistical mechanics, our problem has been 
mapped to an effective Euclidean field theory in d+l dimensions. But, unlike the equilibrium 
case, the parameter t has the meaning of time, not inverse temperature. 

A simple scaling analysis of (||) yields the applicability conditions of the mean field theory 
approximation in d < 2: The rescaling 

t^tT, X ^ y/m^, <f>* (t)^ Zcf), q^Wq (6) 

with ZW = transforms the functional integral measure as follows: 

g-s[0,0%Ao,t.r>] ^ g-^s[0,0M,i,i]^ 

where 

= w 

is the dimensionless coupling. 

We see that the mean field theory (or, equivalently, the weak coupling {g{t) « 1) 
approximation ) is always violated in the limit t — > oo if d < 2. Dimension 2 is critical in a 
sense that the dimensionless coupling g \d=2 doesn't depend on time. 

In the weak coupling limit the average mass distribution n{q, t) solves the Euler-Lagrange 
equation describing the extrema of the effective action (^ : 



dtn{q,t) ^ Y I dq'n{q- q' ,t)n{q\t) - Xq j dq'n{q,t)n{q' ,t), (9) 



in which we recognize the Smoluchowski equation. 

Now we can explain physical meaning of the condition g{t) << 1 of applicability of 
the Smoluchowski equation. Integrating (^) w. r. t. q one finds that the total density 
n{t) = J dqn{q, t) obeys the standard rate equation of the theory of A + A ^ A reaction: 

dtn = -^n\ n(0) = no. (10) 

The asymptotic solution of this equation for t » Uq^Xq^ is n{t) '--^ Note that n{t)~^/'^ 
is the average distance between particles which determines the correlation length and y/lji 
is the typical size of spatial fluctuations of local density. Thus the condition g{t) « 1 is 
equivalent to 



(11) 

l-diff 
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The conclusion is that if the correlation length is much smaller than the scale of fluctuations, 
the spatial inhomogeneity of local density is irrelevant and the mean field theory is applicable. 
Condition ( pi] ) should be compared to the Ginsburg-Landau criterium of applicability of 
mean field theory approximation to the Ising model, see |9| for a review. Note that in 
d < 2 the correlation length measured in units of a typical fiuctuation size diverges as 
i — !■ oo. Spatial correlations become important and mean field theory breaks down in 
complete analogy with the break up of mean field theory approximation of equilibrium 
statistical systems near the critical point. 

Analyzing the applicability of the mean field theory we concentrated on parameters Aq 
and t only, keeping the rest of the parameters of the theory fixed. As a result, (111]) gives 
only a necessary condition of the applicability of Smoluchowski equation. In particular, we 
will see below that the Smoluchowski theory of aggregation can produce incorrect results 
for very small and very large g's even if (|ll] ) is satisfied. 

Smoluchowski equation (H) with the initial condition (||) is easy to solve using Fourier 
transform in q variable with the following result for the average mass distribution: 

"^'^'^^ MF^{t) ^l-F-^i) ' ^ ^ 

where F{t) = 1 + uqIXq. Note that no — > oo limit of (^2|) exists and corresponds to the 
unique self-similar solution of (|^): 

1 g_ 

limna^oon{q,t) = j^^j^e "^o* (13) 

To compute the average mass distribution taking into account the fluctuation effects we 
must analyze the full effective field theory (||) . This turns out to be a relatively simple task 
due to the connection between (|^) and the effective field theory oi A + A which we will 
now describe. 

It is well known that the evolution of total density of particles in constant kernel stochas- 
tic aggregation is governed hy A + A A process. Here we establish a more general result: 
the average mass distribution in constant kernel aggregation is given by the Fourier trans- 
form of the (analytical continuation of) average density of particles in A + A ^ A process 
with respect to the initial average density. 

To verify the statement we note by analogy with jl^ that the last term in the effective 
Hamiltonian (^ is an x- integral of the complete square of the field W = J dqtjffj)- The 
appropriate term in the functional integral measure in (H) can thus be rewritten using the 
Hubbard-Stratonovich transformation: 

g-A / dxdtw^ ^ f £)^(a;^<)e-^ / dxdt{e±^Vxwo_ (^4) 



Substituting (]IJ) into (^ we find that the exponent of the new functional integration 
measure is linear in 0*. The quadratic part of the new effective action is non-degenerate. 
Thus the (/)*-field can be integrated out, at least perturbatively. Following the inverse of 
Martin-Siggia-Rose (MSR) procedure |]ll| one finds that the computation of the average 
mass distribution is equivalent to averaging the solution of the following stochastic integro- 
differential equation with respect to noise: 



Ao 
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{dt-DA)cl>{^,t,q) = ^ dq'(j>{^,t,q-q')cj>{x,t,q') 
-Ao / dq'(j){yi,t,q)(l)(ji,t,q') -2iy/\) ■ Cix,t)(l){x,t,q), (15) 
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where the stochastic potential £,{yi,t) is Gaussian and "white", (^(x, t)^(x', i')) — S{x ~ 
x')(5(t — t'). Imaginary unit multiplying the noise term of ( p^ ) is due to the fact that 
aggregating particles are anticorrelated. Note also that that even though n{q, t) = t)), 
4> and the local mass distribution are different as random variables. For instance, all higher 
moments of (^-distribution are different from higher moments of n-distributon. See pO) for 
more details. 

Rescaling maps (|l5|) to an equation with g{t) entering only as a multiplier of the 
noise term. Thus in the weak coupling limit ( |l5| ) reduces to as it should. 

Equation (^5|) simplifies upon taking a Fourier transform with respect to the mass vari- 
able q: 

{dt-DA)cl){x,t,k) 

= 27rAo0^(x,t, fc) +47rAo0(x,t, fc) + 2i^/)^C{'^,t)(f){-K,t,k). (16) 

Hence the random field 0(x, t, 0) which corresponds to the total density does indeed 
satisfy the closed stochastic partial differential equation (SPDE) arising in the theory of 
A-l-A A reaction, |l^. In general, the field (/)(x, t, k) is coupled only to itself and 0(x, i, 0). 
Introducing new random variables A(x, t) — (/)(x, t, 0) and i?(x, t) — (/)(x, t, 0) — 0(x, t, k) we 
find that they satisfy the following SPDE's: 



{dt - DA)A{x, t) + 27rAoA2(x, t) - 2i VA^$(x, i)A(x, t) = Ao5{t) (17) 
{dt - DA)B{y:, t) - 2tt\qB^{x, t) - 2i^/\'oi{x, t)B{x, t) = BQ5{t), (18) 

where = = f^(l — (1 + ^^)^^) ^re the initial values of random variables A and 

B correspondingly. 

But if g > 0, n{q,t) = - / dke'^^ (B{x,t)) , where B is the solution of SPDE (|l|) of the 
A + A A theory with the initial density parametrised by k. In the universal limit no — > oo 
the initial "density" is 

Ba = ikM, (19) 

thus confirming the statement about the connection between the constant kernel aggregation 
and the theory oi A + A —^ A made earlier. As no flows to infinity under the renormalisation 
group transformation we can set it equal to infinity from the very beginning without 
affecting the answers for i >> Hence everywhere below we assume that Bq is given 



by (19) 



The random field A turns out to be irrelevant for the computation of the average mass 
distribution. Yet it is not completel y d ecoupled from the B-field due to the common ran- 
dom potentials of SPDEs (17) and (|l8|). This remark is important for the computation of 



correlation functions, ||7[. 

The Fourier transform of the average mass distribution n{k, t) is most conveniently com- 
puted within the field-thcorctic formalism. Thus we apply the MSR procedure to ( [l8| ) to 
arrive at the following representation of n{k,t): 



n{k,t)^n{t)- j DBBe^'^^^^+I d^xB'Bo ^ ^20) 
where n{t) is the average density of particles and 

S[B] ^ J^dtJ d'^x{B*{dt - DA)B + 2ttXoB*B^ + Xo{B*B)^} (21) 
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is an effective action of the A + A ^ A theory, Q . 

Now we shall obtain a representation of n{k^t) as a perturbative series in g{t), identify 
the terms in the series giving the leading order contributions in the large t-limit and sum 
these contributions exactly using the formalism of renormalization group. In our analysis 
we are going to use fairly standard methods of the statistical field theory, see e. g. 
for a review. 

The Fourier transform of the average mass distribution can be presented graphically as a 
sum of all Feynmann diagrams of the theory (|2^) with a single outgoing line. These can be 
classified by a number of loops. To regularise the divergencies due to integration along loops 
we use the dimensional regularization, i. e. compute Feynman integrals associated with 
each diagram in dimension d — 2 — where e is a small parameter. A simple combinatorial 
and dimensional analysis shows that a graph with n loops and m ingoing density lines yields 
a contribution to the average density of the form 

C„,™(e)A"-l(i?t)-'^/2-(n-l)(rf-2)/2(^^^^^m^ 

where Cn,m{() is a constant. Therefore, 

oo 

n{k,t) = n(t)^—y^g{trU{e,K), (22) 

where n = —iXqBq = X^iMk is a scaling variable: For a fixed n expansion ( ^ ) is indeed 
an expansion in powers of g{t). An expansion coefficient /„(e, n) is given by the sum of all 
Feynman diagrams containing precisely n-loops. As k is not a small parameter of the theory, 
these sums have to be computed exactly. Direct evaluation of tree and 1-loop contributions 
gives: 



I f m 2/1 Y'^ , , 

Jo (TTT^^ V (2 + e) ~ (2 + e)(4 + e) 

Note that the sum of all tree diagrams yields the solution of the rate equation (|l^) with 
the initial condition (pS|). Hence the expansion ( p2| ) is the standard loop expansion around 
a nontrivial classical solution. 

The inverse Fourier transform of ( p3|) produces the g(t)-expansion of the average mass 
density around the solution (^3|) of the Smoluchowski equation. The coefficients of this 
expansion depend on the mass scaling variable Q = q/{XotM). 

We see that if g{t) << 1, there are domains in n and Q spaces in which the mean field 
theory answer holds. If however d < 2 they shrink with time due to the the growth of the 
dimensionless coupling g{t). Note that we do not discuss the situation in d > 2, where our 
considerations must be modified to account for the strong cut-off dependence of the answers. 

The growth of g{t) is in turn due to the recurrent property of random walks in d < 2: 
the probability that two particles coalesce at some late time t is strongly modified by the 
fact that they met many times in the past and yet didn't react, [p^ . In other words, 
the bubble diagrams renormalising the coupling Aq are all relevant at large times, A 
simple dimensional argument shows that only the renormalisation of two particles processes 
is relevant at large times in d < 2. Thus one can extract the large time asymptotics of the 
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average mass distribution from (^2|) by rearranging the 3(t)-expansion in such a way that it 
becomes an expansion w. r. t. a renormahsed dimensional coupUng gnit). This is precisely 
the task the renormalisation group method has been designed for. 

The contributions corresponding to diagrams renormalising Ao form a geometric sequence 
and thus can be summed exactly. The result is 

Mt) = r, (24) 

where Cd is a positive constant. Suppose that the average mass distribution at i = 
is known. Then we can in principle compute it for any later time t by solving the field 
theory (|2l|), but with Aq and no(fc) replaced with Ai{(to) and n{k,to) correspondingly. The 
result for t > to will not depend on to, as the evolution operators U{t',t") corresponding to 
( pT] ) form a semi-group: C/(0,io) • U{tQ,t) = U{0,t). With the aid of dimensional analysis 
condition dt„n{k^t) = can be written in the form of Callan-Symancyk (renormalization 
group) equation: 



(25) 



Here 



f3{g{t)) ^ -t^^ = CMtf eg{t) (26) 
is an exact /3-function of the field theory (^l|). The solution of Calan-Symancyk equation is 

n(fc, t)=(^^y (^^^ ^ k,to, gRit, to)) , (27) 
where g{t, to) is the running coupling constant: 

gM{l + ^lni±)y' : d = 2 (27.2) 



Expressions (|2^) and (g7|) can be used in the analysis of the large time asymptotics of 
the average mass distribution in the following way: 

If d < 2, limt_>oo 5i?Xt, to) — which is a non-trivial fixed point of the /3-function (p^). 
Thus the large time asymptotics of the average mass distribution in c? = 1 has according to 
( p7| ) the following self-similar form: 

n{q,t) = ^<^>( (28) 

where <& is some universal scaling function. Hence our analysis reproduces the exact scaling 
properties of the average mass distribution in one dimension known from the exact solution, 
[ p^ . The substitution of (^) into ( p7| ) gives a representation of n(g, t) in the form of e-series. 
Unfortunately, the e-expansion doesn't reveal any information about the scaling function $ 
when e = 1 (d = 1). 

If however d = 2, then gR{t,to) ^ An/lnit/to) « 1 for t » to- Hence it is possible 
to evaluate the r. h. s. of (p7|) for large times using the weak coupling expansion (p3|)! 



7 



Technically, this possibility is realized as follows: We express g{to) in the series (^3|) through 
fffl(io) using ( p^ ) for e > and substitute it into (p7|). The resulting expression is non- 
singular at e = 0, as the renormalization of the coupling constant Xq takes care of all 
perturbative divergencies of the field theory (|2^). Thus we can set e = and after taking 
the inverse Fourier transform obtain the following answer for the average mass distribution 
in d = 2: 

(2(1 - QR)iln{QR + C) + 2Qn + l- ^{Qr - 2)2)) + 0(.gj,(i, <o)')), (29) 

where C is the Euler constant and Qn — MgH(t to)-Dt renormalized scaling variable. 

The cutoff time can be related to the lattice spacing: to = Ax^ /D. Even though individual 
terms in the expansion ( p9| ) depend on to, their sum doesn't. A simple check shows that 
integrating (p9[ ) with respect to q gives the one-loop approximation of the total density 
obtained in |p|7while the integral of q-n{q, t) gives the total mass density M. The expression 
for the average mass distribution can be simplified even further if we use a more natural 
scaling variable Q = q ■ n(t)/M, where n{t) is the exact total density. Then the term 
S3^Aol(2C — 5) (2 — Qfl) disappears from the r. h. s. of ( |29| ) and we are left with the 
non-trivial form corrections to the scaling function e^^. 

Thus, we conclude that if Qr » to/t and Qr « {ln{t/toY^'^ , 

^ M\ 4:7: Dt J ^\ inDtM J' ^ ' 

thus establishing the result obtained numerically in jj] and conjectured on the basis of 
heuristic arguments in 1^ . Note that (|3^) holds at large times irrespectively of the value of 
the coupling Xo/D. But if the latter is large, the mean field theory result ( psf ) isn't valid at 
any time at all. 

We emphasize the non-uniformity of the corrections to the renormalized mean field an- 



swer (30) with respect to the scaling variable Qr. This phenomenon occurs in higher dimen- 
sions as well. A computation shows that in d < 4, the large-Q asymptotics of the one-loop 
correction to the mean field result behaves like g{t)Q^e^'^ . Thus the mean field theory can 
be violated even above the critical dimension if Q is large enough. 

The main reason for the violation of (|3^) for large g's is the existence of another length 
scale in our problem, Im = {q/M)^/'^ (so far our analysis of the applicability of the mean 
field theory used only the diffusion scale Id — Vot). Corrections to the mean field theory 
which depend on Im are due to the fluctuations in the flux of masses in the mass space due 
to aggregation, see for more details. 

In conclusion I would like to stress that the applicability of the statistical field theory 
methods to the stochastic aggregation is not necessarily limited to the simple model at hand. 
In particular, it would be interesting to apply the machinery of renormalization group to 
the model of stochastic aggregation with a multiplicative kernel K{q,q') = qq' , for which 
the mean field theory predicts a finite time gelation transition, [p^ . Another potentially 
interesting application is to a model of stochastic aggregation with fragmentation which 
admits a non-equilibrium phase transition between a trivial and a non-trivial steady states, 
[|6| 
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